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Abstract: The measurement of sparticle masses in the Minimal Supersymmetric 
Standard Model at the LHC is analysed, in the scenario where the lightest neutralino, 
the x?j decays into three quarks. Such decays, occurring through the baryon- number 
violating coupling X"- k , pose a severe challenge to the capability of the LHC detectors 
since the final state has no missing energy signature and a high jet multiplicity. We 
focus on the case A 212 7^ which is the most difficult experimentally. The proposed 



method is valid over a wide range of SUGRA parameter space with A 2 



212 



10' 



0.1. 



Simulations are performed of the ATLAS detector at the Large Hadron Collider. 
Using the Xi from the decay chain qi — > X2Q. ~^ Ir^Q ~^ Xi^Qi we show that the 
Xi and X2 ma sses can be measured by 3-jet and 3-jet + lepton pair invariant mass 
combinations. At the SUGRA point m = 100 GeV, m 1/2 = 300 GeV, A = 300 GeV, 
tan/3 = 10, fj, > and with A 2 ' 12 = 0.005, we achieve statistical (systematic) errors 
of 3 (3), 3 (3), 0.3 (4) and 5 (12) GeV respectively for the masses of the X11 
and qL, with an integrated luminosity of 30 fb _1 . 
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1. Introduction 

The addition of supersymmetry (SUSY) to the Standard Model (SM) represents a 
theoretically attractive way of addressing several of the problems faced when at- 
tempting to reconcile constraints from fundamental models at high scales with the 
phenomenology seen at the electroweak scale. In general, supersymmetric models 
which attempt to solve the naturalness problem, and fix the Higgs mass near the 
electroweak scale, predict a rich spectrum of physics which can be accessed by exper- 
iments at the Large Hadron Collider (LHC). However, most studies of experimental 
signals for SUSY have assumed that R-parity (Rp) is conserved. Rp is a multiplica- 
tive quantum number defined as (— 1) 3B+L+2S where B and L are baryon and lepton 
numbers, and S is the spin of the particle. It therefore distinguishes SM particles 
(Rp = +1) from their superpartners (Rp = — 1). 

If Rp is conserved (RPC models), SUSY particles can only be created in pairs, 
and the lightest SUSY particle (LSP) is stable. Therefore SUSY events each contain 
an even number of LSPs, which escape detection and give rise to a large missing 
transverse energy (E™ 1SS ). This signature has been exploited by many analyses of 
the discovery potential of the LHC [1], since it provides a clean separation between 
SUSY events and the SM background. However, the incomplete measurement of the 
final state makes the reconstruction of the SUSY mass spectrum more difficult. 
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The following terms [2] may be added to the Minimal SUSY Standard Model 
(MSSM) superpotential in order to incorporate R-parity violation (RPV): 

Wft, = ^XijkLiLjEk + X' ijk LiQjD k + -X'^UiDjD,, + mUH^ (1.1) 

where gauge indices have been suppressed. 

In the MSSM Lagrangian, SM and SUSY particles are grouped together into 
lepton (Li), quark (Qj) and Higgs (H 12 ) SU(2) doublet superfields and electron 
(Ei), down (Di) and up (f/j) SU(2) singlet superfields. In Equation |L~j] , Xijk, A- jfc 
and X'lj k are Yukawa couplings between the matter superfields and /tj is the mixing 
term between the lepton and Higgs doublets. The subscripts i, j and k are family 
indices. X^k is antisymmetric under i <-> j and A" jfc is antisymmetric under j <-> k. 
Wft p therefore adds 9 + 27 + 9 + 3 = 48 free parameters to the MSSM superpotential. 

Interactions involving the A^, A- jfc and Ki couplings violate lepton number, while 
those involving the X"j k coupling violate baryon number. The simultaneous presence 
of the second and third terms in Equation [LI] can lead to fast proton decay in gross 
conflict with the lower limit on the proton lifetime [3] . Since both lepton and baryon 
number violation are required in order for the proton to decay, current experimental 
bounds on the proton lifetime and other SM processes can be respected if either 
baryon number or lepton number is conserved [2]. 

The difference in experimental signatures between RPV and RPC SUSY models 
at the LHC depends on the strength of the RPV coupling. We will concentrate on 
the trilinear couplings and neglect the bilinear term which leads to mixing between 
the leptons and gauginos, and between the sleptons and Higgs bosons. When the 
RPV couplings are small compared to the MSSM gauge couplings, the dominant 
effect is that the LSP can decay into SM particles. For example, the lifetime of the 
LSP (x?) as a function of the RPV coupling A 2 'i2 is shown in Figure [T|a. 

If the RPV couplings and MSSM gauge couplings are of the same order of mag- 
nitude, RPV production processes and decays of particles heavier than the LSP 
become important. For a large X'^ k coupling, the branching fractions of RPC and 
RPV decays of a squark can be of the same order of magnitude (Figure |l]b), and 
SUSY particles can be produced singly, as was investigated in [4]. 

In the present analysis, X% 12 is the only RPV coupling set to a non-zero value. 
This gives rise to the LSP decay mode Xi ~~ > °ds. This is the most challenging case 
experimentally, since there are no leptons or 6-quarks among the xi decay products 
which can be used as tags for signal events. As each event will usually contain two 
Xx' s there will be at least six jets in the final state. 

The RPV coupling was added to a minimal supergravity (mSUGRA) model with 
5 GUT-scale parameters: a universal scalar mass mo = 100 GeV, a universal gaugino 
mass 7711/2 = 300 GeV, trilinear Hff soft SUSY breaking terms Aq = 300 GeV, the 
ratio of the vacuum expectation values of the two Higgs doublets tan /? = 10 and the 
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(a) Lifetime 



(b) Branching Ratio 




Figure 1: (a) Lifetimes of the d,R and Xi an( f (b) branching ratio of RPC (dashed) and 
RPV (solid) decays of the cLr, plotted against A^ at the SUGRA point too = 100 GeV, 
m 1/2 = 300 GeV, A = 300 GeV, tan/3 = 10 and sgn/x +. 

sign of the SUSY higgsino mass parameter \i positive. It should be noted that we 
have only included the RPV coupling at the weak scale, i.e. the RPV coupling is not 
used in the evolution from the GUT scale, as was done in [5]. 

Five sets of parameter values have been extensively studied in the Rp conserving 
MSSM. The parameters chosen here correspond to SUGRA Point 5, with one mod- 
ification: the value of tan/3 has been increased from 2.1 to 10 in order to keep the 
predicted Higgs mass above the current experimental limit. The masses of some key 
particles in this model are given in Table [|. Searches in the SUGRA Point 5 scenario 
have been well studied in the case of a stable Xi [!]• At this SUGRA point, the xi is 
the LSP, as must be the case for our analysis to be valid, even though cosmological 
constraints which require the LSP to be neutral only apply if it is stable. 

The case of A^ = 0.005 is considered first. This coupling strength gives rise to 
decay chains essentially identical to those in an RPC model, except for the decay 
of the xi inside the beam pipe, with a lifetime of 1.0 x 10~ 14 s. Unlike many other 
RPV couplings, A^ is n °t currently constrained by experiment [6]. 
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Table 1: Masses of selected particles (GeV) for the model investigated. 



The effect of varying the RPV coupling in our analysis will be discussed in 
Section || For RPV couplings of order 10~ 6 or smaller [7], the LSP has a sufficiently 
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long lifetime to decay outside of the detector. If the LSP is charged, heavily ionizing 
low velocity tracks would be seen in the detector, providing a clear signature. This 
analysis addresses the case of a neutral LSP, the lightest neutralino (xi), which has 
negligible interactions with the detector. For low RPV couplings, the experimental 
signature is then identical to that of an RPC model. However, if RPV couplings 
are above 1CT 4 [7], the LSP usually decays in the beam pipe and the missing energy 
signature, seen in RPC models, is not present. 

2. Analysis Strategy 

In this work, HERWIG 6.1 [8] is used as the event generator 1 , and the official ATLAS 
simulation program ATLFAST [10] is used to simulate the performance of the ATLAS 
detector. Since each xi decays to three quarks, and in general the decay chain 
produces the xi i n association with at least one other quark (typically from squark 
decay), the mean jet multiplicity (A/jet) in the signal events is high. 

The principal difficulty in measuring the xi mass is the identification of the 
correct jets from the xi decay. Nearly all right-squarks decay via — > XiQ ~^ QQQQ 
and one might therefore expect A/j ct = 8 for <jRqR production. Gluon radiation by 
quarks, however, raises this to an average of 9.2 jets, in spite of the fact that the three 
jets from harder x°s are spatially close together and some merging of jets occurs. In 
QlQl events Aj et = 10.7. The increase with respect to the right-handed states is due 
to the difference in couplings to charginos and neutralinos. Gluinos mostly decay 
into a squark and a quark and gg events have a higher value of N- iCt = 12.8. A simple 
algorithm is used to calculate the jet energies, summing the energy within a cone of 
size 0.4 about the jet axis in the r\ — <f> plane, and at least 8 jets with Et > 25 GeV 
are required in signal events. 

The analysis proceeds in the following steps: 

• Cuts are applied to reduce the SM background to below 10% of the SUSY 
signal. These cuts rely on the presence of lepton pairs in the signal events. 
Such lepton pairs are produced from the decay chain x° — > Zr Z — > llxl i n most 
SUSY models. An analysis based on looking for events with two such X2 decays 
was proposed in [11] but the rate for events with four leptons is much lower 
than for events with only one x\ decay of this type. 

• In each signal event, cuts are made on the jet transverse momenta (pt) in order 
to preferentially select jets from neutralino decays. 

• All possible combinations of three of the selected jets are inspected, and their 
invariant mass, rrijjj, calculated. Events are retained if two combinations are 
compatible with the same candidate xi mass. 

lr The simulation of RPV events in HERWIG is discussed in [9]. 



4 



• One of these three-jet Xi candidates is combined with an opposite sign, same 
flavour (OSSF) lepton pair. The invariant mass of this system, rrijjju, is a xl 
candidate. A clear peak at the X11X2 masses in the rrijjj, rrijjju plane is then 
observed. 

• The l R and the masses are reconstructed using 3-jet combinations in the 
2-dimensional Xi, X2 mass peak. 

In the following sections, each of these steps is considered in turn. 



3. Standard Model Background 

The SM background in this model is considered in [1]. It is shown in [12] that the 
inclusive SUSY signal can be separated from the SM background by requiring that 
each event contains: 

• at least 8 jets with Et > 25 GeV; 

• at least one jet with Et > 100 GeV; 

• transverse sphericity> 0.2, transverse thrust< 0.9; 

• ^T,cent > 1 TeV, where mr,cent = Yl Pr* + S P^^i where the sum includes 
central jets and leptons (i.e. with pseudorapidity \q\ < 2); 

• at least two leptons (e or //) with p T > 15 GeV and \rj\ < 2.5. 

With these cuts, the signal to background ratio is greater than 10. The SM 
background has not been explicitly simulated in this study. Current Monte Carlo 
event generators are not capable of reliably simulating QCD eight jet plus two lepton 
production. We have therefore simulated eight jets and two leptons distributed 
according to phase space. 

In SUSY events, lepton pairs are created in the decay X2 W — * HXv The 
leptons are therefore required to have opposite charges and the same flavour (OSSF). 
The invariant-mass distribution of the lepton pairs created in this decay has a kine- 
matic edge [1, 13] which is given by 

m „ / K(x°) - m\l R )] x [m\l R ) - m\x^\ ,„ n 

mee =\ , (3.1) 

V m 2 (l R ) 

and is simulated after experimental resolution in Figure |2|. 

With the particular parameter set chosen this edge is calculated as 95.1 GeV. 
Accordingly, events are required to have a lepton pair with an invariant mass below 
this value. The corresponding edge can be easily measured for other parameter sets. 
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Figure 2: The dilepton invariant mass for events with an OSSF electron or muon pair, 
after the SM cuts have been applied. The kinematic limit for the decay chain shown in 
t igure | is at 95.1 GeV. Events are excluded if there is no jet combination which passes the 
jet cuts described in Section |I[ For the SUGRA point chosen it happens that the kinematic 
edge lies just above the peak at m(Z°) from the decay Z° — ► li. 



4. Detection of the Xi an d xi 

Many different decay chains can contribute to the SUSY signal selected by the cuts. 
One important example, qi — > X2Q. ~^ ~^ Xi^Q ~ > QQQMQi is shown in Figure [| 




Figure 3: One of the decay chains of the <jl contributing to the signal. 

When A212 is small, there are nearly always two x? s produced in an event and 
one can therefore search for two sets of three jets with similar invariant mass. An 
upper limit on the invariant mass difference of Srrijjj = — rnf--\ < 20 GeV is 

used in this analysis, where a and b label the two Xi LSP candidates. 

In order to limit the combinatorial background, the search for the Xi signal is 
initially restricted to events with 8 < A/j et < 10, with the following cuts on the 
allowed range of jet transverse momenta in GeV: 
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10 20 30 40 50 60 70 80 90 100 
jet p T / GeV 



Figure 4: The rescaling factor applied to the observed jet energy as a function of jet pt- 

• 100 < z4 al) ; 17.5 < pf ] < 300; 15.0 < ^ a3) < 150; 

• 17.5 < # } < 300; 17.5 < pf 1 < 150; 15.0 < pf 1 < 75, 

where al denotes the highest pr jet from neutralino candidate a, etc. 

Candidate sets of jets from the xi decay can also be identified by their separation 
in the rj — plane. For both Xi^, cuts are made on the distance between the hardest 
and next hardest jets (AR12) and between the combined momentum vector of the 
two hardest quarks and the softest quark (A_R 12 _3). The following cuts are chosen 
based on simulations: 

• AR$ < 1.3; A/^-s < 1-3; 

• AR^ < 2.0. 

Since the SUSY cross section is dominated by production of squarks and gluinos, 
about 95% of events have two hard jets with E^ l) > 200 GeV and E!f 2) > 100 GeV 
from the squark decays. We require that two jets in an event satisfy these cuts. We 
do not use those two jets to construct Xi candidates. This decreases the background 
from wrong combinations. 

For each combination of jets passing the kinematic cuts, the jet energies are 
rescaled according to their pr, to allow for energy lost out of the jet cones. The 
rescaling function, shown in Figure |J is the ratio p^/Pt used in ATLFAST. 

The reconstructed masses of all Xi candidates are included in Figure ^. If two 
or more combinations of jets from one event pass the cuts, the masses from each 
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m... / GeV m... / GeV 

J J J tliJtP 

(a) (b) 

Figure 5: (a) The invariant mass of three-jet combinations passing the cuts described 
in the text. The mass peak from the decay x? — * Q11 can be seen above the background 
from wrong combinations of jets. The input Xi m ass is indicated by the arrow, (b) A 
phase-space sample shows a peak in much the same region. 



combination are plotted with unit weight. However this can lead to one event con- 
tributing a large number of combinations and therefore if there are more than five 
combinations which pass the cuts only the five with the smallest difference between 
the Xi candidate masses are included. 

In events with more than one combination passing the cuts, two combinations 
of jets often differ only in the choice of jets for one of the two Xi s - I n such cases, the 
unique mass combination is included in the histogram only once. The two ambiguous 
masses are both included with unit weight. 

There is a broad combinatorial background beneath the x\ mass peak in Fig- 
ure |5|a, the shape of which is defined by the kinematic cuts and is reproduced by the 
phase-space sample, as shown in Figure |5|b. In order to further suppress the back- 
ground we attempt to find the mass of the x® m the decay chain X2 — > Ir I ~ y MXi 
by forming the total invariant mass of the OSSF dilepton pair and one of the 3-jet 
candidates. We choose the Xi candidate which is nearest in rj — to either lepton. 

By using the extra information from the leptons we are able to supress the 
combinatorial background. A clear peak in the (X11X2) m &ss plane is visible in 
Figure |6|a. Figure [F| shows slices through the peak along the axes m ± = {rrijjj ± 
m jjju) I V%- The peak is present in the central slices, while the sidebands show similar 
shapes to the background under the peak. It is clear that this peak is not determined 
by the kinematic cuts, as it is absent in our phase-space sample (Figure B|b). In 
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Figure 6: (a) The Xi { m jjj) an d X2 ( m jjju) candidates. The number of jet combinations 
passing the cuts per 30 fb _1 is given by the key. The circle and the ellipse show the peak 
and standard deviation of a 2-d gaussian fitted to the data contained in the dashed box. 
The star shows the input masses, (b) The corresponding (iTijjjjTnjjju) invariant mass 
combinations from the phase-space sample show no such peak. 



addition we find that the position of the peak accurately follows the input masses, 
when those are varied from (116.7, 211.9) to (137.8, 252.6) GeV. 

The data in the rectangle shown in Figure |6] were fitted with a 2-d gaussian. Since 
rrijjj and rrijjju are highly correlated, the peak was fitted in the rotated (m + , m~~) 
coordinate system in which the correlations are smaller. The mass difference relies 
on lepton rather than jet momenta, so the width in the m~ direction (4 GeV) is 
smaller than in the m + direction (15 GeV). The standard error on the peak was 4.5 
GeV in m~ and 1.5 GeV in m + . This corresponds to a 3 GeV uncertainty in each of 
the neutralino masses. 

Our fitted masses, at m(xi,X2) — (118. 9, 215.5) GeV are slightly high when 
compared to the input values of (116.7, 211.9). This is due to several effects, in- 
cluding overlap between the jets, and the contribution of energy from the underlying 
event in the jet cones. Indeed we would expect that our simple rescaling factor will 
overcompensate for energy losses from jet cones, since the three jets from the Xi 
decay are close in 77 — 0, so energy losses from one cone can end up in one of the 
other two. 

A fuller investigation of these effects is beyond the scope of this study, requiring 
as it does a full investigation of the calorimeter calibration procedure for multi-jet 
events. It is estimated that in the actual experiment, the uncertainty in the absolute 
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Figure 7: Slices taken through Figure in the rotated (m + , m~) coordinate directions. 
The fit was performed on the data contained in the region between the dotted lines. The 
fitted peak and width are given by the thick line and arrow at the top of each plot. 

jet energy scale will be of the order of 1% for jets with pt > 50 GeV [1]. For lower 
energy jets an uncertainty of 2-3% is more likely. With real data, therefore, it may 
be possible to reduce the systematic uncertainty on the Xi mass to 3 GeV. 

5. Detection of the Ir and qi 

For the measurement of the slepton mass, having fitted the Xi an d X2 masses, we 
select combinations within 1 x o of the peak. These combinations, with two OSSF 
leptons preferentially select the decay chains x\ ~^ ltd MXu an d X2 ~^ Z°Xi- Only 
the former can be used for Ir measurements. The region of mSUGRA parameter 
space in which this chain will exist is given by the condition m(xl>) > m {^R)i an d is 
shown in Figure || 

The dilepton invariant mass distribution, Figure ^, for our point in MSSM pa- 
rameter space shows only a very small peak at the Z° mass, but a clear kinematic 
edge, indicating that the slepton decay chain x_2 IrI — > Hxl dominates. This is 
expected at this point, since the m(%2) — m (Xi) mass difference of 95.2 GeV means 
that there is little phase space available for the decay X2 ~ * Z°X?- If the masses and 
couplings for this decay were significant, then we would exclude events with dilepton 
invariant mass near m(Z°) from the In measurement. Both chains may have been 
preceded by qi — > X2Q1 so both can be used for the qi mass measurement. 

We recall that the initial sample of Xi candidates was restricted to events with 
8 < iVjct < 10 in order to reduce the combinatoric background. However, choosing 
combinations from under the (X11X2) mass peak removes much of the background, so 
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Figure 8: The region of mSUGRA parameter space over which the decay x° — * Ir£ ~~ * Xi^ 
occurs is shaded. Mq is the universal scalar mass, and M1/2 is the universal gaugino mass 
at the GUT scale. The lower hatched region is excluded by lack of electroweak symmetry 
breaking or the existence of tachyonic particles, whereas the in the upper hatched region 
the Xi is n °t the LSP and our analysis does not apply. The contours show the mass of the 
ul squark. The x marker shows the chosen point. The other parameters are: tan/3 = 10, 
A = 300 GeV and sgn(^) = +. 

in this section the jet multiplicity cut is therefore relaxed to 8 < A/j ct < 11 in order 
to increase the statistics. The invariant mass of the three-jet neutralino candidates 
is adjusted to the best-fit mass of the x?, by rescaling the Xi jet momenta by the 
same factor. 

The Ir mass is found by combining the Xi candidate closest to a lepton in 
77 — <p with that lepton. The resulting invariant mass distribution, rrijjji, is shown 
in Figure |^a. The sharp peak was fitted with a gaussian, with another gaussian for 
the background. This gave m(ln) = 157.8 ± 0.3 GeV, which is slightly high when 
compared to the input value of 155.3 GeV for the same reasons as were discussed in 
Section |j. 

The experimental electron and muon momentum scale uncertainties are expected 
to be small fractions of 1% [1], so the systematic error in the slepton mass measure- 
ment will be dominated by the same (3 GeV) jet scale uncertainty as m(x 1 ). The 
statistical error in rescaling the 3-jet invariant mass to the fitted m(xi) peak in- 
troduces another 3 GeV systematic error into the Ir and q~R masses. The overall 
systematic error in m(7#) is therefore 3 © 3 = 4.2 GeV. 

The hardest two jets in the event are assumed to have come from squark or 
gluino decay, and so are not used in making three-jet combinations for neutralino 
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Figure 9: The masses of (a) Ir and (b) <jl candidates. Only combinations from under 
the (X11X2) m ass peak, and satisfying the cuts described in the text are plotted. Each jjj 
invariant mass has been rescaled to the fitted Xi mass. The Ir and qL masses are indicated 
by the arrows. The fitted functions are described in the text. 

candidates. The mass is found by combining each x% candidate with the harder 
of these two leading jets. To reduce the background we select combinations within 
2xcT of the m(Z#) peak. This allows us to relax the cut about the (X11X2) P ea k from 
1 to 2 x a. The resultant invariant mass distribution, rrijjjuj, is shown in Figure |9|b. 
A peak is visible near the ul and dt masses of 633 GeV and 638 GeV respectively, 
but the resolution is not sufficient to separate the states. 

The background was modelled by finding the invariant-mass distribution of the 
X2 candidates with the hardest jet from other q L candidate events. A gaussian fit to 
the signal with this background shape gave m(qi) — 637 ± 5 GeV. The uncertainty 
in modelling the background was estimated by fitting the distribution with other, 
simpler background shapes. These decrease position of the peak to 634 GeV (for a 
resonance-shaped background) to 627 GeV (for a linear background). This shows a 
systematic uncertainty in the qL mass of about 10 GeV. The hard jet used in the 
calculation of m(qi), haspr > 100 GeV introducing an uncertainty in the mass scale 
of 1% [1], or 6 GeV. Carrying forward a 3 GeV uncertainty in the jjj invariant mass 
scale and 3 GeV from the Xi fi% the total systematic error in the squark mass is 
12 GeV. 

At this SUGRA point the dominant decay mode of the q~R is qR — > XiQ- One 
might therefore try to reconstruct the q~R mass by combining the Xi candidate not 
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used in the reconstruction with the second hard jet. However while the cuts we 
applied to reduce the Standard Model background, i.e. requiring the presence of 
two leptons, mean that all the signal events contain a (Jl, they do not necessarily 
contain a right squark. Only a third of the signal events actually contain a q~R. Most 
of the SUSY events at this SUGRA point come from gluino production which will 
either be rejected due to the large number of jets or contain additional hard jets 
from the gluino decay, and hence have a large combinatoric background for the (Jr 
reconstruction. In an attempt to reduce this background it is possible to use cuts on 
the X2 an d Ir masses from q~i decay on the other side of the event such that there 
is only one c[r candidate. However this reduces the statistics so that a signal cannot 
be observed. This combination of factors makes it impossible to reconstruct the q~R 
mass at this SUGRA point with the available statistics. 



6. Other values of 

The method outlined above is relatively insensitive to the size of the coupling A^- 
However as the RPV coupling A^ is decreased, the lifetime of the Xi increases as 
shown in Figure |l|a. The method will start to fail when x^s predominately decay 
beyond the first tracking layer of the detector. Special reconstruction could increase 
this by about an order of magnitude, at which point the decays would occur outside 
of the tracking volume. We therefore exclude events when one or other of the Xi s has 
travelled more than 100 mm (1000 mm) from the interaction point in the transverse 
direction. 

As can be seen in Figure ffU] statistics become limiting for A 2 ' 12 less than about 



10~ 5 , when cr « 800 mm. With smaller couplings the RPV decay of the Xi effectively 
switches off, and a RPC analysis based on a missing transverse energy + lepton(s) 
signature, such as [13], is effective. 

If A212 is larger than 0.1, an initial q~R often decays immediately into 2 jets 
and then only one 3 jet invariant mass combination will necessarily be close to the 
neutralino mass. 

The size of the RPV coupling can be determined from the Xi lifetime, as shown 
in Figure |I]a. The lifetime is, in principle, measurable for a wide range of couplings, 
by using vertexing information. However the need for detector-level Monte-Carlo 
simulation makes the measurement of A^ beyond the scope of this paper. 



7. Conclusions 

The usual ubiquitous signature for the Rp conserving MSSM is missing transverse 
energy. This signature disappears once Rp violating couplings are added, unless they 
are extremely small and the lifetime of a neutralino LSP is such that it decays outside 
the detector. 
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Figure 10: The number of jet combinations with transverse decay lengths less than 
100 mm (thick line), and less than 1000 mm (thin line), as a function of the RPV coupling. 
Combinations beyond 1 x a of the (XiiXlD peak found for A^ = 0.005 are excluded. 

We examined the case where the neutralino LSP is unstable and decays to 3 jets 
with no particular tags on them (for example b's). This corresponds to the RPV 
coupling X'212, the trilinear RPV coupling giving the hardest case in which to detect 
and measure sparticles. 

By analysing the decay chain, qi — > X2Q ~^ Ir^Q — > Xi^Qj we have shown that 
the Xii X2 an d Ql be detected and their masses measured, and that the mass of the 
Ir can also be obtained in much of parameter space. The sparticle production and 
decays in this signal are all Rp conserving apart from the x? decay into 3 jets. 

Although we have used a point near SUGRA Point 5 to derive the soft SUSY 
breaking parameters, the method should in principle work for other more general 
SUSY breaking parameter sets in which the decay chain in Figure ^ exists. When 
some of the sparticles involved in the decay chain become much heavier than 1 TeV, 
the analysis will become statistics limited. 

To summarise, in the MSSM with a trilinear RPV coupling, even in the hardest 
choice of A^, it is possible to detect sparticles and measure their masses at the LHC. 
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